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Reheating or entropy production on the three-brane from a decaying bulk scalar field is studied in the brane-world picture 
of the Universe. It is shown that a significant amount of dark radiation is produced in this process unless only the modes 
which satisfy a specific relation are excited, so that subsequent entropy production within the brane is required in general 
before primordial nucleosynthesis. 
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Any new theory of gravity and/or high energy physics must pass a number of cosmological tests, among which is im- 
plementation of cosmological inflation pip). Successful inflation must fulfill three requirements, namely, sufficiently long 
quasi-exponential expansion driven by vacuum-like energy density such as a potential energy of a scalar field, termination of 
accelerated expansion associated with entropy production or reheating to set the initial state of the classical hot Big Bang 
cosmology well before the primordial nucleosynthesis and generation of primordial fluctuations with desired amplitude 
and spectrum 0. It is much more difficult to achieve the second element than the first in general, and the third one typically 
requires fine tuning of model parameters. 

In this paper we consider reheating after inflation in the brane world picture of the Universe [^[^). In this scenario, our 
Universe is described on the four-dimensional boundary (three-brane) of ./^-symmetric five-dimensional spacetime with a 
negative cosmological constant A5 = —6k 2 , where k is a positive constant. This situation not only takes into account the 
spirit of Horava-Witten theory J^Us), but also recovers the Einstein gravity around the brane with positive tension |^PJiq[. 
Much work has been done on brane-world cosmology fll|-|l3| including inflationary brane solutions Q-^. 

We assume the five-dimensional Einstein gravity with a negative cosmological constant A5 and a three-brane at the fifth 
coordinate w — about which the spacetime is Z2 symmetric. We write the metric near the brane in the following form in 
terms of the Gaussian normal coordinate. 

ds 2 = gABdx A dx B — —N 2 (t, us)dt 2 + Q 2 (t, w)a 2 (t) (dx 2 + dy 2 + dz 2 ) + dw 2 = q^dx^dx" + dw 2 , (1) 

where capital Latin indices run 0,1,2,3, and 5 while Greek indices from to 3. We take N = Q — 1 on the brane w = 0. One 
can write down functional forms of N(t,w) and Q(t,w) explicitly as 



2 (t,w) = cosh(2A:w) + \k~ 2 H 2 (cosh(2fc™) - 1) - ^1 + k~ 2 H 2 + Ca~ 4 sinh(2fe|w|), 



N\t,w) = Q~\t, 
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in the case the bulk is in a vacuum state with a negative cosmological constant A 5 , where C is an integration constant J13|]. In 
this solution the induced metric on the brane is nothing but the spatially flat Robertson- Walker metric with the scale factor 
a(t). 

The evolution equation on the brane in this case is given by 

A 4 k\ 2 k 2 C 



2 ( a\ K50- A4 K 5 2 k C , . 
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A 4 . - ^A 5 + -fo' j , (4) 

where k 2 is the five dimensional gravitational constant related with the five dimensional reduced Planck scale, Ms, by 
k 2 — A/g -3 . Here a is the brane tension, p to t is the total energy density on the brane, and the last term of (^) represents 
the so-called dark radiation with C being an integration constant 0,0jl3). In order to recover the standard Friedmann 
equation with a vanishing cosmological constant at low energy scales, we require a = 6k/ k, 2 and k 2 = k|ct/6 = ft 2 fc, where 
k 2 is the four dimensional gravitational constant related with the four dimensional reduced Planck scale, M4, as k 2 — M^ 2 . 
We therefore find Mf = Mf/fc. That is, if we take k = Ma, all the fundamental scales in the theory take the same value, 
k = M4 = Ms. Note that k also sets the scale above which the nonstandard term quadratic in p to t is effective in (^|). We 
assume that k is much larger than the scale of inflation so that such quadratic corrections are negligible. 

We consider the case inflation is driven by a bulk scalar field (j> with a five-dimensional potential V[(j>] and study the 

evolution of <f> after brane inflation, because reheating is expected to proceed in the same way as in four dimensional theory 
if the inflaton lives only on the brane . Since <j> is homogenized in three space as a result of inflation, it depends only on t 
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and w. We consider a situation that 4> rapidly oscillates around = and assume that V[4>] is expressed as V[4>] = m 2 (j} 2 /2. 
Then the Klein-Gordon equation reads 

□ 5 <M<M) - V[<j>(w,t)] = -j=d t {^g ao 4>) + -^=d w (V=^,„) - V[0] = 0, (5) 

where a dot denotes time differentiation. In order to express energy release of (j> we introduce the following dissipation term 
phenomenologically in (^|). 

a^( w ,t)-V'mw,t)] = ^H^ + r B I^. (6) 

Here To and Tb represent energy release to the brane and to the entire space, respectively. The denominator in the right- 
hand-side is introduced on dimensional grounds. 

From (^) and (^) together with the Z2 symmetry, we find 

£i» = -£i = ^(0,*). ( 7 ) 

where superscripts + and — imply values at w — > +0 and —0, respectively. 
The divergence of the energy-momentum tensor of the scalar field, 

T% = <t>,u<t>,N - g MN (h PQ ^p^Q + V[4>]) , (8) 

reads, 

T { % = W(»,t) - V'\4>(w,t)]} <P, A 
Integrating A — component of ([]) from w = — e to w = +e near the brane, we find (Q) from the zeroth order in e and 
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with 



-{3H + r B )4> J (o,t)-Mo,t), (10) 



^ = _^ 2 + V[0], ^ = -^^(7=^/), (11) 
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from the terms proportional to e. Thus the energy dissipated by the Fd term on the brane is entirely compensated by the 
energy flow onto the brane. 

Next we study how the energy released from (j> affects evolution of our brane Universe by analyzing gravitational field 
equations [ fui|]l7| ] . In the present situation the total energy momentum tensor including the contribution of bulk cosmological 
constant reads 

Tmn = —K,5 2 A 5 g M N + T^flj + SmnS(w), (12) 

where Smn is the stress tensor on the brane. Its nonvanishing components can be further decomposed as 

Sfiv = — aq^v + Tfiv. (13) 

Here t^ v represents energy momentum tensor of the radiation fields produced by decay of <j> and it is of the form t£ = 
diag(— p r ,p r ,p r ,p r ) with p r = p r /3. 

In terms of the unit vector tim = (0, 0, 0, 0, 1) normal to the brane, the extrinsic curvature of ate =constant hypersurface is 
given by Kmn = lM ( l < N n Q;P with Iain = <?a/jv — Umtin- Then from the Codazzi equation and the five dimensional Einstein 
equation, we find 

D V K^ — D^K = n 2 5 T M Nn N = n^T^ = K5<j><j>, w 5°, (14) 
where D v stands for the four dimensional covariant derivative with respect to the metric q^ v . The above equation reads 

D V K" + -D K+ = 4^^(0,1), (15) 

near the brane, w — > +0. 

From the junction condition and the .^-symmetry, on the other hand, we find 

K+, = -4(^-l<l^), (16) 
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therefore 



D V K V + - D,K+ = -^D.S; = -^-D vT ;. (17) 



Combining ([L5|) and (|l^), we obtain 



r 



DutZ = ~4>% t (18) 



namely, 



^ = SH(Pr +Pr) + |f 2 = -4^ + (19) 



on the brane. Thus we find only the dissipation term proportional to To with the delta function is effective to reheat the 
brane. This equation has the same form as the reheating in perturbation theory after conventional inflation in four dimensional 
theory 

On the other hand, the four dimensional Einstein tensor, G^t", satisfies the following equality on the brane Jl7| , 

G ( f = k\ (T«f + r;) + 4< - E%, (20) 

with 



T 



(*)- 
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(21) 



Here 7r^ represents terms quadratic in r„ which are higher order in p r l{kM^) 2 and are consistently neglected in our analysis. 
E^ = C*™^ is a component of the five dimensional Wyel tensor Cp ! q , which is the origin of the dark radiation jl3[ . 
Now we write down the four dimensional Bianchi identity, 

D u &*> =Q = kI (D v T^ )v + D v tZ) - D u E;, (22) 

to yield 

2 



D V K = - 



K 4 d fl - 2 1 ,2 . 1 2,2\ 2k 4 H - 2 K 4 r )2 / OQ \ 

2fc^ " 2^ + 2 m ^ J ~ ^ - 2^ ' (23) 



Note that <f> 2 w = Y 2 D 4> 2 / '(16k 2 ) on the brane and is negligibly small compared with <p 2 , since we expect the dissipation rate 
To is at most comparable to the scale of inflation and hence much smaller than the scale k. With this approximation the 
quantity in the parenthesis in ( pj|) may be replaced by pj,(0,t). We also find D v Eq = BqEq + 4HEq because is traceless. 
As a result we obtain the following set of evolution equations in the brane universe w = in terms of tp(t) = c/>(0,i)/y/2k. 

a\ 2 «4 / , . \ _ 1 .2 . 1 2 2 P<p _ E o , nA \ 

-J = y (P^ + Pr + e), P V =-^P +2 m <P = 2fe> £= "^f' ( ^ 

^ = -(3if + r fl ^ 2 -J^ J,^, (25) 
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4ii-p r + r D ^, (26) 



= -4ffe-(if + r D -r B )^ + J ¥ „ (27) 



where we have used (px|). 

Here J^, represents energy flow of the scalar field into the bulk and we cannot calculate it unless we solve the field equation 
in the bulk. It is, however, formidable to do so in the realistic situation that the bulk metric is influenced by the bulk scalar 
field which would be the case in brane inflation driven by a bulk scalar field. In this case the metric would be even more 
complicated than (Q). Hence in order to obtain some insights on the form of J v we must employ several approximations and 
find a bulk scalar solution. 

Let us assume the bulk metric is governed by A 5 and further neglect terms suppressed by fc _1 and Ca~ A in (|j|). Then the 
metric reads 

ds 2 = - e ~ 2k ^dt 2 + e - 2fcW a 2 (i) (dx 2 + dy 2 + dz 2 ) + dw 2 , (28) 

and field equations (||) and (|^) become more tractable. In the time scale of the field oscillation the dissipation terms are 
unimportant because they are perturbatively small quantities by assumption. So let us estimate J v by solving (^) instead of 
(|) first. The solution of (§) under the metric ( Bjj ) is easily found in the same way as Goldberger and Wise Jl^] as 
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with 



T n (t)^a-^(t)e- im "\ 

Y n (w) = e 2ftM [j„ (^ e fc l-l) + bn N v (^e fcH )] 



v = 2 



1 + 4k 2 ~ 2 + 4k 2 : 



(30) 
(31) 



under the assumption that the field oscillates rapidly in cosmic expansion time scale. Here m n is a separation constant which 
may take continuous values in the present case with a single brane, and b„ is a constant determined by the boundary condition, 
= 0, at w = as 



^(t) + ^(t)]Mt) + t*(t 



(32) 



If we incorporate the effect of dissipation on the boundary condition, 



it is modified to 



2 + 



im n r D 
2k 2 



J " ( T ) + m [( 2 + * (t) + T« (x) 



(33) 



where use has been made of T n (t) = —im n T„(t). 

Hereafter let us assume only a single oscillation mode exists. This is a natural assumption because any higher mode is 
expected to decay earlier and the final stage of reheating would be dominated by the lowest mode that has been excited. 
Then we find 



(m n - m jtptp. 



(34) 



Since oscillates rapidly in the expansion time scale let us average the right-hand-side of evolution equations (|2jj)-(|27|) 
over an oscillation period. Using <p 2 (t) = m\ip 2 (t) and <p<p(t) = —(3H + YB)<fi 2 (t)/2, we obtain the following set of evolution 
equations in the brane universe w = where a bar denotes average over the oscillation period. 

2 „2 



(!) =y(^+^ + £ ). 
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(35) 
(36) 
(37) 
(38) 

(39) 



We assume that both p r and e are vanishingly small at the end of inflation when we set the initial condition, because they 
rapidly redshift during inflation and very little creation is expected in that period |l7| . We identify the reheating epoch with 
the time when p v becomes smaller than p r . From ( pi^ ) this occurs at t ~ H^ 1 ~ '" 
place during H >Yb- Here we qualitatively analyze the system for various values of m 



r B , so that the dominant creation takes 



Case A: m„ > m 

In this case, the creation terms of dark radiation, — (Yd + H — rejin^ 2 + ^(3H + YB)('rn 2 — m^ip 2 are negative and 
their magnitude is larger than the creation term of radiation, Yovn^ip 2 during the important period H > Yb- Hence 
more dark radiation is created than ordinary radiation in magnitude. Then p r and e tend to cancel each other and 
the higher-order terms of the Friedmann equation, which we have neglected so far, would play an important role in the 
subsequent evolution of the three brane. This means that we do not recover standard cosmology on the brane after 
inflation. 

Case B: m n <C m 

This possibility is not excluded in the brane-world scenario with non-factorizable geometry unlike in the case of ordinary 
Kalza- Klein compactification with factorizable metric In this case the last term of (S) is dominant and we find 

more dark radiation than ordinary radiation unless Ts is extremely small with Yb/Yd < m n /m 2 <C 1. 
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Case C: m n ^ m 

This is the most delicate case and the final amount of dark radiation can be either positive or negative depending on 
the details of the model parameters. But in order to achieve successful primordial nucleosynthesis, the amount of extra 
radiation-like matter is severely constrained Ji). In order to have sufficiently small e compared with p r after reheating 
without resorting to subsequent entropy production within the brane, the magnitude of creation terms of e should be 
vanishingly small at the reheating epoch H ~ Ts. That is, only the mode that satisfies the inequality 

|2r s (m 2 -ml) -T D m 2 n \ < r D m 2 n , (40) 

should be present. 

We thus find a specific relation (^) should be satisfied for the graceful exit of brane inflation driven by a bulk scalar 
field <j>- In case it is not satisfied, we must introduce some other mechanisms of entropy production within the brane before 
primordial nucleosynthesis which imposes stringent constraints on the exotic energy density and the expansion law at that 
era One way is to assume that cj> predominantly decays into massive particles whose energy density redshifts less rapidly 
and dominates over e soon. If these particles decay into radiation before nucleosynthesis creating an appropriate amount of 
baryon asymmetry, we may recover the standard cosmology. Another way is to consider second inflation which is driven by 
a scalar field confined on the brane |l^| to dilute e. 

In summary we have studied entropy production on the three-brane from a decaying bulk scalar field <j> by introducing 
dissipation terms to its equation of motion phenomenologically. We have shown that the so-called dark radiation is significantly 
produced at the same time unless the inequality (jlo|) is satisfied. Although we have analyzed only the case with a specific 
form of the dissipation, we expect our conclusion is generic and applicable to other forms of dissipation, too, because it is 
essentially an outcome of the four dimensional Bianchi identity (E2I). We therefore conclude that in the brane- world picture 
of the Universe, it is likely that the dominant part of the entropy we observe today originates within the brane rather than 
in the bulk. 
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